Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 February 2008 (MN M£X style file vl.4) 



On the contribution of density perturbations and 
gravitational waves to the lower order multipoles of the 
Cosmic Microwave Background Radiation 

A. Dimitropoulos^ 1 )* and L.P.Grishchuk^ 1 ' 2 ) t 

(1) Cardiff University, Department of Physics and Astronomy, P.O.Box 913, Cardiff, CF24 3YB, United Kingdom 

(2) Sternberg Astronomical Institute, Moscow University, Moscow 119899, Russia 

7 February 2008 

ABSTRACT 

The important studies of Peebles, and Bond and Efstathiou have led to the formula 
Ci — const/ [£(£+1)] aimed at describing the lower order multipoles of the CMBR tem- 
perature variations caused by density perturbations with the flat spectrum. Clearly, 
this formula requires amendments, as it predicts an infinitely large monopole Co, and 
a dipole moment C\ only 6/2 times larger than the quadrupole C%, both predictions in 
conflict with observations. We restore the terms omitted in the course of the derivation 
of this formula, and arrive at a new expression. According to the corrected formula, 
the monopole moment is finite and small, while the dipole moment is sensitive to 
short-wavelength perturbations, and numerically much larger than the quadrupole, 
as one would expect on physical grounds. At the same time, the function £(£ + l)Ce 
deviates from a horizontal line and grows with £, for £ > 2. We show that the inclusion 
of the modulating (transfer) function terminates the growth and forms the first peak, 
recently observed. We fit the theoretical curves to the position and height of the first 
peak, as well as to the observed dipole, varying three parameters: red-shift at decou- 
pling, red-shift at matter-radiation equality, and slope of the primordial spectrum. It 
appears that there is always a deficit, as compared with the COBE observations, at 
small multipoles, I ~ 10. We demonstrate that a reasonable and theoretically expected 
amount of gravitational waves bridges this gap at small multipoles, leaving the other 
fits as good as before. We show that the observationally acceptable models permit 
somewhat 'blue' primordial spectra. This allows one to avoid the infra-red divergence 
of cosmological perturbations, which is otherwise present. 
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1 INTRODUCTION 



The important studies of Peebles (Peebles 1982; Peebles 1993) and Bond and Efstathiou (Bond & Efstathiou 1987; Efstathiou 
199C) have led to the formula for the multipole moments of the CMB radiation, Ce, caused by small density perturbations 
with flat spectrum (also known as Harrison-Zel'dovich-Peebles, or scale-invariant spectrum; spectral index n = 1): 

const 



Ce 



(1) 



This formula aims to describe the lower order multipoles of the CMBR variations in a spatially flat universe. The general 
acceptance of this result is reflected in the fact that the Ce distributions are usually plotted in terms of the function £(£+ l)Ce 
versus £. In these units, formula (|l|) is conveniently plotted as a horizontal line. 
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Clearly, formula (|l|) was derived under some assumptions, and this derivation requires amendments. Indeed, let us start 
from the monopole moment Co, setting £ — in this formula. Expression (|l|) predicts infinity for Co, in conflict with 
observations. Surely, it would be difficult to observationally distinguish a monopole component caused by perturbations from 
the isotropic unperturbed temperature To, but it is strange that a theoretical formula gives infinity for an effect produced 
by small perturbations. Let us now consider the dipole, £ = 1, and the quadrupole, £ = 2, moments. According to (|l|), the 
quantities Ci and Ci are numerically of the same order of magnitude, Ci = (6/2)C2. However, we know from observations 
that the observed dipole variation of the temperature is about 100 times larger than the observed quadrupole variation. 

The predictions of formula (llj) also seem to disagree with independent theoretical studies. According to Grishchuk & 



Zel'dovich ( 1978), in the limit of long- wavelength perturbations, which normally provide dominant contributions to the £ = 
and 1 = 2 multipoles, the monopole and quadrupole variations of the CMBR temperature should be small and approximately 
equal. At the same time, the dipole variation should be suppressed, as compared with the monopole and quadrupole, in 
proportion to the ratio la/X, where In is the Hubble radius and A is the perturbation's wavelength, A ^> lu- The multipole 
moments Ci, being quadratic in the temperature variations, contain this factor in the second power. The factor £h/A becomes 
much larger than 1 in the opposite limit of short-wavelength perturbations, A <C In - The dipole component is, thus, expected to 
be very sensitive to short waves, and its value should eventually exceed the values of the monopole and quadrupole components, 
as soon as the contribution of short-wavelength perturbations is also included. However, formula (|l|) suggests results differing 
from these conclusions. The discrepancy is especially troublesome, as the derivation of (^) seems to have taken into account 
the perturbations of all wavelengths. Indeed, this formula arises from the table integral (discussed below in Section 

V I t\l2 dn 1 

— 



where the integration over wave-numbers n ranges formally from zero to infinity. Therefore, according to this integral, the 
very short- wavelength perturbations have also been taken into account in the calculation of all C^'s, including the dipole Ci. 

The demonstrated deficiencies of formula (|l]) for £ = and 1 = 1 indicate that there could also be significant deviations 
from the values suggested by (|l|) in the case of £ = 2 and higher order multipoles. This consideration has prompted us to 
reanalyse the derivation of formula (|l[). In Section U we introduce the necessary notations and remind the reader the basics of 



the Sachs-Wolfe result (Sachs & Wolfe 1967). For the purpose of later comparison with observations, we consider primordial 
spectra with arbitrary spectral indices, so that the flat spectrum n = 1 is simply a particular case. Section |^ summarises 
the statistical assumptions and defines the multipole moments. In Section ^ we identify the terms that were dropped out 
in the course of the original derivation of formula (jij). We show the importance of these terms and restore them. After this 
correction, the numerical values of the monopole and quadrupole moments become comparable to each other, while the dipole 
moment becomes sensitive to short-wavelength perturbations, as it should. At the same time, the correct formula for Ci 
requires the function £(£ + l)Ct to deviate from the horizontal line and grow with £ for £ > 2. Formally, this growth would 
be unlimited if one were allowed to use the primordial power-law spectrum for arbitrarily large values of the wavenumber n. 
However, it is known that the primordial power spectrum of density perturbations should experience a turnover, caused by the 



transition of the cosmological expansion from the radiation-dominated to the matter-dominated regime ( Zel'dovich & Novikov 



198J 



Peebles 1993). By the time of decoupling, the primordial spectrum gets modified. Mathematically, the original power 
spectrum gets multiplied by the modulating function (transfer function). The modulating function is practically unity in the 
spectral region of long-wavelength perturbations, but it strongly suppresses the power in short-wavelength perturbations. As 
a result, the growth of the function £(£ + l)Ce terminates and is being followed by its decline at multipoles, roughly, around 



£ = 200. This causes the function £(£ + l)Ce to form the first peak, recently observed (Hanany et al. 2000 
et al 



2000 



de Bernardis 



Mauskopf et al. 200C). The following (less prominent) peaks of the modulating function are responsible for the 
further peaks of the function £{£ + l)Ci. In Section ^ we analyse in detail the form of the modulating function and derive 
the possible (theoretical) positions and heights of the first peak. We think that the origin of the subsequent peaks in the 
£-space is due to the modulation of the metric power spectrum of density perturbations (oscillations in the n-space). This 
modulation, in its turn, entirely arises due to the standing-wave character of the initial density perturbations. Travelling 
sound waves in the radiation-dominated era would not have produced Sakharov oscillations at the time of decoupling, as was 



explained long ago (see Section 11.6 in Zel'dovich & Novikov (1983)). The recent paper (Miller et. al. 2001) may be providing 



support to this interpretation of the peaks. However, the origin and position of secondary peaks is not the subject of this 
paper, so we leave the door open for alternative interpretations. The full computation of all multipoles is presented in Section 
^| We replace the existing uncertainties about the matter content of the Universe by the two phenomenological parameters: 
the decoupling red-shift z^ cc and the matter-radiation equality red-shift z eq . We also vary the third parameter, (3, which 
characterises the spectral slope of the primordial perturbations. (The flat spectrum n = 1 corresponds to j3 = —2.) We show 
how the distributions of multipole moments depend on the choice of these parameters. Since we concentrate on anisotropics 
of gravitational origin, the chemical composition of matter in the Universe is unimportant, as long as the perturbations can 
be described by solutions to the perturbed Einstein equations at the radiation-dominated stage (equation of state p = e/3) 
followed by the matter-dominated stage (equation of state p = 0). Section (7] gives a preliminary comparison of the calculated 
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multipole moments (including the dipole) with the available CMBR observations. We require the theoretical (statistical) 
dipole, as well as the location and the height of the first peak, to be as close as possible to their observed values. It appears, 
then, that the values of the COBE-observed (Smoot et al. 1992) multipoles, around £ = 10, are somewhat larger than the 
theoretical curves for £(£ + l)Ce would require. The, so far neglected, intrinsic anisotropics (temperature variations at the 
last scattering surface) can only increase the height of the first peak making this discrepancy larger. The most natural way of 
removing this discrepancy is to take into account gravitational waves. They mostly contribute to the multipoles from £ = 2 to 
approximately £ = 30, without changing the dipole, and not seriously affecting the location and height of the peak. Section 
|§] describes the contribution of gravitational waves. In effect, we calculate the amount of gravitational waves needed in order 
to raise the 'plateau' at lower £'s to the actually observed level. In Section ^ we again compare the theory with observations, 
this time including gravitational waves, and discuss the values of the parameters Zdcc, z cq , f3 that produce distributions 
£(£ + l)Ce in better agreement with CMBR observations. We argue that a slightly larger amount of gravitational waves allows 
the primordial spectrum to be somewhat 'blue' (n > 1 or ft > —2) instead of being flat. The CVs produced by such 'blue' 
spectra are still consistent with the discussed observations. We consider this fact as a satisfactory conclusion, since, from the 
theoretical standpoint, the flat spectrum and 'red' spectra (n < 1 or f3 < —2) are undesirable. Indeed, the mean square values 
of the metric (gravitational field) fluctuations are logarithmically divergent for the n = 1 case and power-law divergent for all 
n < 1, in the limit of long waves. The removal of these divergences would require extra assumptions. Appendix A contains 
some technical details of the derivation of the corrected formula for Ce- 



2 THE SACHS- WOLFE RESULT AND GROWING DENSITY PERTURBATIONS 



Following Landau & Lifshitz (1975), Sachs & Wolfe (1967), we describe a FLRW universe and the perturbed gravitational 
field by the line element 

ds 2 = a 2 [dr/ 2 — (Sij + hij)Ax l dar'], 

where hij refers to density perturbations or gravitational waves. In the matter-dominated era (equation of state p = 0) the 
scale factor behaves as a(r]) oc rf , which we write as 

a m (rj) = 2Z H (v - ?7m) 2 , 

where In is the Hubble radius at the present time r\ = 7/r. The constant r) m is needed for further considerations. It is convenient 
to choose tjb. — T]m = 1. With this convention, the wave-length equal to the Hubble radius today corresponds to the wave- 
number riH = 4-7T. In the matter-dominated era, the employed coordinate system is both synchronous and comoving - the most 



convenient choice (Landau & Lifshitz 1975, Sachs & Wolfe 1967). The elements of plasma emitting the CMB photons, as well 



as the observer, are 'frozen' in the perturbed deforming matter and, by the definition of the comoving coordinate system, 
have zero velocities with respect to the matter fluid. 

Let the direction of observation of a photon emitted at 77 = t/b be characterised by a unit vector 



e = (sin# cos<f>, sin9 sin<f>, cos6 



(2) 



The relative temperature variation 8T/T seen in a given direction (9, <f>) is described by the integral (Sachs & Wolfe 1967) 



ST 1 
T (e) = ^ 



drj 



dw, 



(3) 



where £ = t^r — tjb and the integration is performed along the path r\ = r)R — w, x = etu. The quantity £ is related to the 
red-shift z^ cc at decoupling as 



(i + z dc .y 1/2 . 



(4) 



Typically, the quantity £ is very close to 1. The fact that the photons coming from different directions might have been 
emitted at slightly different tie's is irrelevant for this calculation, as the integrand in expression (^) is already a small quantity 
of first order. A first order variation of the upper limit of integration will change the result only in second order of smallness. 
However, a small initial variation of the temperature over the last scattering surface will be reflected in an additive term to the 
integral (^). This term would be effective even in the absence of the intervening gravitational field hij along the photon's path, 
and it can be called the intrinsic anisotropy (8T/T)i n - We discuss this term later, in Section [| As long as further scattering 
and absorption of CMB photons can be neglected, the integral (^) plus the (8, 0)-dependent initial condition (5T/T)in fully 
describe the temperature variations at the point of reception. 

We now consider density perturbations. By the time of decoupling, one can normally neglect one of the two linearly- 
independent solutions to the perturbed equations, and work with the so-called growing solution. In the case of density 
perturbations at the matter-dominated stage, the growing solution is 
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tain,*) = b(x)^- + iv ~^ m)2 m,^. 

Using formula (jE|) in the integral (|^), one obtains (Sachs & Wolfe 1967) 
—(e) = — B f e f \ - — (1-0 B f e f \ + —B\ R - —B\ E . 



10 



'B 10 



(5) 



(6) 



10 1 10 

The exceptional property of the solution (Q) is that the integrand becomes a total derivative, so the integral (j3|) reduces to 
the values of the under-integral quantity at the end-points of integration. This property is responsible for formula In the 
general case, this is not true, and is only approximately true for other types of perturbations when they are taken in the 
long- wavelength limit (Grishchuk & Zel'dovich 1978). In view of the already existing disorder in the literature with respect 
to the naming and interpretation of different parts of the Sachs- Wolfe work, we simply call the terms in (g) , in order of their 
appearance, as term I, term II, term III, and term IV. Term I is purely dipolar in its angular structure. Term II consists of 
all spherical harmonics, including the monopole. Term III has no angular dependence: it only contributes to the monopole. 
Term IV consists of all spherical harmonics, including the monopole. 

One normally assumes that cosmological perturbations can be decomposed over spatial Fourier components with arbitrary 
wave- vectors n. In our case, this amounts to the decomposition of -B(x) as 



B(x) 



(2tt) 3 / 2 



d3 ( Ti in-x . j-i* — in 
n ( B n e 4- B n e 



'— 2 J d 3 n B n {in ■ e [l - (1 - Z)e^] + (l - e~-*) } 



Since -B(x) is a real function, we have included a complex conjugate part to make the expression manifestly real. Using the 
notation 

- (2^') 

the temperature variation (0) can be written as 
f<e) = ^[G{e)+G*(e)]. 

One can expand G(e) over complex spherical harmonics: 

oo I 

G(e)=J2 12 9imY em (e,<p). 

1=0 m=~l 

To do this explicitly, it is convenient to write the wave vector n in the form 
n=(nsinO cos$, nsinO sin«I>, ncosO) 



(8) 



(9) 



(10) 



and use the expansion (see, for example, formula (16.127) of Jackson ( 1975)): 

oo i 

exp{in • x} = 4tt ^ 12 ^ ^ nX ^ Y «™( 6 ' *) Y ^{6, 4>)- 



(11) 



t=0 m=-t 

As a result, one finds that 
term I produces a purely dipolar variation 



1 oo £ 

in . e = i -£ n ^ ^(O, *) Yij($, cb) = 4tt ^ ^ i e \n 6n F/ m (6, *) Y im (9, 



term II contains all spherical harmonics 



in • e (1 - e in e « 

1=0 m=-l f=—l 

term III is a purely monopolar contribution 



^ 12 12 12 Y n(1 ~^ Jf(<) Yi *f(®'^ r ^( e ' $ ) Y ^(°> 



1 = 4vr Y w {6, <t>) Y *<£Q, *) = 4tt J2 12 ' Sm y ^ e ' $ ) Yim (°> < 

1 = m=-l 

term IV contains all spherical harmonics 

oo I 

exp{in -eO=4*1212 ' «("0 Y ^ 8 ' ^ Y ^ 9 ' 
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The part of gi m produced by the sum of terms I, III and IV is easy to calculate: 



(i) 



4tt 



(2tt) 3 / 2 



d 3 n B n 



The part of gi m contributed by term II demands a more tedious calculation, which results in 
(2) 4tt . 



(2tt) 3 / 2 



d 3 n B n n(l - f) 1 



21 + 1 



[(i + 1) - £ ^_iK)] y/ m (e, <&). 



An outline of the calculation may be found in Appendix [a|. The total expression gi„ 
4ir 



(i) (2) 

= 9lm+ 9lm IS 



d J n B n D,(n) F/ m (9,$), 



(2tt) 3 / 2 

where De(n) denotes the combination 

D t {n) = in S a + n(l - 1^ + L ) Jt+W ~ 4 J<-i("01 + 5 « - 



(12) 



(13) 



The terms in Di(n) shall be referred to as term A, term B, term C and term D, in the order of their appearance in (|13[). They 
originate from the terms I - IV introduced earlier, respectively. For the three cases I — 0, £ = 1 and I > 2, De(n) is explicitly 
given by: 



A>(n) = n(l-C)jiK) + l-io«), 

Di(n) = in + n(l - £) ± [2j 2 «) - jo«)] - 



(14) 



^> 2 (n) = n(l - 



2/:- 



-[(£ + 1) j f+1 «) - t jt-frQ] ~ Hni). 



Finally, we can write 

00 1 

^ = ^ Y t Jfi, 4>) + 9tm YfJO, 0)] (15) 

1=0 m=-t 

where the ge m are given by ([L?]), 

The calculation of the integral requires only the knowledge of the gravitational field perturbations hij. However, for 
further references, we remind the expression for the growing density contrast (5 pj p)( x i V) associated with the solution (|B|): 

The density contrast is a product of a purely n-dependent part (77 — rj m ) 2 and a purely spatial part 



(16) 



Equation ( |16[ ) and the decomposition 

" -f- OO 



f (x) = (2*)-"" 



d n 



+ I — 

P 



establish the link between the Fourier components of the density variation and the gravitational field perturbation: 



(17) 



(18) 



3 STATISTICAL ASSUMPTIONS AND MULTIPOLE MOMENTS 

Formulae (^) , (^) give a temperature variation over the sky, which is produced by a single and deterministic, even if very com- 
plicated, perturbed configuration of matter density and gravitational field. We believe, however, that primordial cosmological 
perturbations may have had a quantum-mechanical origin. To this end, the quantities hij and 8T/T are quantum-mechanical 
operators, and the extraction of the observable information should proceed through the calculation of quantum-mechanical ex- 
pectation values. Later on, we shall go into some details of the quantum-mechanical generation of cosmological perturbations. 
However, in this Section, we want to keep the discussion at a phenomenological and as simple as possible level. 

We make a simplifying assumption that the quantities B n are classical random functions. Since B n appears in the 
expression of 5T/T, the temperature variation is also a random function. The angular brackets < ... > will denote the ensemble 
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average. A single statistical hypothesis that we shall need in the subsequent calculations is expressed by the equalities (B n ) — 
and 

(B n B;> = |B n | 2 S\n - p) = (B p B*), (B n B p > = = (B* n B*). (19) 
The quantity \B n \ 2 characterises the power spectrum of -B(x). Indeed, by multiplying (^) with itself and using (|l^), one finds 

(B(x) B(x)) = / ^n 3 |B„| 2 . (20) 

The quantity we need to calculate is the correlation function 
/ST. , <5T , 

The vectors ei, e 2 refer to two different directions of observation, separated by the angle 712. Using the product of two 
expressions (^), one finds 



(^r(ei) Y"(e 2 )) = YM^I X! ^ 9lm g * l ' m '} Ytm ( 01 ' ( l >1 ) Y £'m'{02,fa) + (dim 9l>m>) Y t Jfi\,<l)i) Y t t m r{$2,<h) 

1=0 m=-(l'=0m'=-J' 

+ {gtm 9 I'm') YlJfil^x) Y e * m ,(9 2 , 4> 2 ) + (gtm gt'm') YlJfi-t^l) Y{/ m /(02, <te)] • (21) 

Using (p^), (p^[), and the orthogonality condition for spherical harmonics (see, for example, formula (3.55) of Jackson ( 1975)), 
one derives 

{gtm gt'm') = = {gl m g\i m i) 
and 

{gem ge>m>) = {gem gt' m ') = 27r _1 {-if i e+e 



(I 



— n 3 |B„| 2 B,(n) D t ,{n) 
n 



Sit' 8mm>- (22) 



The substitution of (|22[) in (|2l[) and the application of the addition theorem for spherical harmonics (see, for example, 



formula (3.62) of Jackson (J 1975)), brings the correlation function to the form 



ST ST sr^ 21 + 1 

< y-(ei) — (e 2 )> = 2^ °e ^(cos7i 2 ), 

where the multipole moments C'e are given by the expression 

Ce = ^rSr f — n " |B "! 2 1) ' (n) - (23) 

100 J n 

Thus, the multipole moments are fully determined by the dimensionless power spectrum n 3 |B„| 2 and the square of the known 
function (|l3|). 

There is no physical reason why the function |B n | 2 should be a power law function of the wavenumber n, let alone to 
have a fixed spectral index in a broad interval of n. For cosmological perturbations of quantum mechanical origin, a power 
law spectrum in terms of n is generated by a power law (in terms of rf) scale factor a,i(rf) describing the very early Universe. 
But, every deviation of Oiirj) from a power law in r\ produces a deviation of the generated spectrum from a power law in 



n (Grishchuk & Solokhin 1991). Virtually any observed |B n | 2 as a function of n can be attributed to a specially chosen 
generating function a\(ri). Had we known |B n | 2 a priori, we would find the other cosmological parameters by comparing Ci 
with observations. If we knew the other cosmological parameters with certainty, we would deduce from observations | i?,T. | 2 
and ai(rf). In reality, however, almost everything should be found from observations. To simplify the problem, we shall follow 
the common tradition and postulate the power law dependence 

|B„| 2 =K{P) n 2 ' 3+1 , (24) 

where the coefficient K{f3) is a constant and j3 is a number. In the theory of quantum mechanical (superadiabatic, parametric) 
generation of cosmological perturbations, f3 enters the power law scale factor Oairj) = lo\ri\ 1+l3 , with j3 < — 1. But for the 
purposes of this presentation, /3 can be regarded as a phenomenological parameter, describing the spectral index. The constant 
K(/3) is calculable from the quantum normalisation, but for the purposes of this presentation it suffices to consider it simply 
as a number. 

It is customary to write the power spectrum of (5p/p)(x) in the form 
P(k) oc k n , 

where k is our wavenumber n. We denote the spectral index by a roman-style n, in contrast to the wavenumber, which is 
always denoted by an italic n. It is easy to relate n with 0. Indeed, using (^), jl9| and (H), one finds 
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(x) 



(x) 



n 3 (n \B n \ 2 ). 



According to our notation, T\k) is (n 4 |_B n | 2 ), so the spectral index n is related to f3 by the relation n = 2/3 + 5. The flat 
spectrum is defined by the requirement that the dimensionless r.m.s. metric perturbation n 3/,2 |i?„| (per logarithmic interval 
of n) does not depend on n. This gives /3 = — 2 or n = 1. 
The substitution of (§3) in (pjh leads to 



Ce = 



4-7T 



100 



■KQ3) 



i 2 ^ +2) Dj{n] 



dn 
n 



(25) 



This formula reduces, for I > 2, to an expression similar to equation (17) of Abbott & Wise (1984) (given there without 
derivation) . 



4 FORMULA (1) AND THE LIMITS OF ITS APPLICABILITY 



Formula (|l|) arises from the (ill-justified) assumption that all terms in (g) can be neglected with respect to the last one, term 
IV. Accordingly, only the last term in (|l3|), term D, is being retained. Formula (|2f]) is, thus, truncated to 

-i 



C e = 



4n 



100 



K{(3) 



-MOV 



dn 

n 



(26) 



We use a calligraphic Ce to denote the multipole moments resulting from this assumption, in contrast to the Ce that follow 
from the correct formula (p5|), where all terms are being retained. Using standard relations (see paragraph 10.1.1 of Handbook 
of Mathematical Functions ( 1985) and formula (6.574.2) of Gradshteyn & Ryzhik (198C)), the integral (Ejj|) is analytically 
calculated ([Peebles 1982} [Bond fc Efstathiou 1987|): 



Ce = 



i 2p+1 ,-208+2) IT-2/3-2) ITI + /3 + 2) 

25 [rx-/3-i)] 2 rx^-/3) w ' 



for ^> -0-2 and /3 < -1. 



For the particular case of the flat spectrum, the multipole moments are found by substituting (3 
(507T)- 1 K{-2) 



-2 in ©: 



c 



for £ > 0. 



(27) 



(28) 



^ + 1) 

This expression is in fact formula (0), the const being (507r) _1 K(— 2). Strictly speaking, the monopole moment ^ = is not 
contained in (pg|), but the suggested infinity follows from a rigorous calculation. One needs to return to (p(i|), setting 1 = 
and j3 = —2, to derive 



C = 



47T" 



100 



K(/3) 



poo 


sin(n^) 




[ «) J 



dn 

n 



(29) 

0, the quantity 



In the short wavelength limit n — > +oo, the integral (g9|) is well behaved. However, in the opposite limit n 
sin(n£) / '(n£) ~ 1, and Co is determined by the logarithmically divergent expression. 

The predicted infinity of the monopole moment is not a reflection of some inherent property of the flat spectrum, or some 
artificial 'gauge problem', as many think, but is strictly the result of using ( Ej| ) instead of the correct formula (Ej|). As was 
already mentioned in the Introduction, formula Jjil ) also predicts that the dipole, £ = 1, and quadrupole, 1 = 2, moments are 
of the same order of magnitude, in conflict with observations. 

It is instructive to explore the integrand of formula (|2f]) separately in the long-wavelength n <C rm and short-wavelength 
n 2> nu regimes. This analysis will show when and what kind of dangers one can expect. It will also demonstrate some 
dramatic differences in the asymptotic behaviour of the integrands in ( |25| ) and (pj^). The most spectacular modifications occur 
in the long- wavelength behaviour of the monopole moment and short- wavelength behaviour of the dipole moment, but all 
multipoles are seriously affected. For this analysis we shall use a strictly power-law spectrum (p^), bearing in mind, though, 
that the transfer function effectively suppresses the power in short waves, so that a real spectrum is turned down at sufficiently 
large n's. It is convenient to separately examine the cases t = 0, I = 1 and I > 2 in each regime. 

The long- wavelength contribution to Co: The expression of D({n) is given by the first line in (|l4|). The participating 
terms are B, C and D. We shall show that the presence of the purely monopolar term C cancels out the infinity arising due 
to term D, making the monopole moment finite and small. Expanding the Bessel functions jo(n£) and ji{n£) up to (n£) 2 and 
(n£) respectively, one derives the approximate expression for D<{n) in the long wavelength limit: 



D (n) 



6 
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Term C, which is 1, exactly cancels with the leading order term in the expansion of term D, which is 
is of the order (n£) 2 <C 1 and not unity. Formula (Eq) yields 



Co 



4ty 



100 



2(0+4) 



-1. As a result, D<{ri) 



(30) 



It is clear that the monopole moment Co, unlike Co, is finite in the case of the fiat (y3 = — 2) spectrum. The integral (|30j) is 
convergent for all /3's satisfying the constraint 

P > -4. (31) 

The monopole moment is genuinely divergent in the long- wavelength limit only for such 'red' spectra, that /3 violates (|jl|). 

The long- wavelength contribution to Ci: D\(n) is given by the second line in (|lj). The participating terms are A, 
B and D. The purely dipolar term A combines with terms B and D in such a way, that the dipole moment in this regime is 
suppressed, as compared with the monopole and quadrupole. With the help of appropriate expansions of the Bessel functions 
one derives 



Din) 



1 



i n£+ -n 3 e 
3 ^ 30 ^ 



(3 - 2^ 2 n 3 



3 3 v 30 v 3 * 30 * 30 

Term A combines with the leading order terms in B and D. As a result, D\(n) is of the or der (n£) 3 and not (n£), as would 
be the case if only term D were retained. In agreement with Grishchuk & Zel'dovich 
(which is a small number for long waves) as compared with Ddn) and D-^n) (for D^n) see (fe 



1978 ), Di(n) has one extra power of n 
below). Formula ([ 



20Y 



dn 



2(/3+5) 



The long- wavelength contribution to C^> 2 : The expression for Dt>^n) is given by the third line in (jwj) 
participating terms are B and D. Expanding the Bessel functions, one finds 



Di>in) 



K+*(i-0]r 



1 • 3 • . . . • (2£ + 1) 
and formula (|25[) yields 



47T 



100 



■K(f3) 



K + ^(i 



1-3- 



dn 



2(0+2+1) 



yields 
(32) 
The 

(33) 
(34) 



As expected ( Grishchuk fc Zel'dovich 1978| ; |Grishchuk 1994 ), Co, given by (^), and C2, given by ([ii]) for I — 2, have exactly 
the same behaviour in the limit n — > 0. The integral (|3j) converges if /3 satisfies the constraint /3 > —2 — I. Not surprisingly, 
this constraint for 1 = 1 reduces to (|3l|). The same constraints on /3 follow from the requirement that the quadrupole moment 
produced by gravitational waves does not diverge in the long- wavelength regime (Grishchuk 1993). 

We now examine formula ( pq 1 ) in the short wavelength limit. In this analysis, we ignore the fact that the primordial power- 
law spectrum is bent down by the modulating (transfer) function and, in any case, the integration should be terminated at 
some large n — n max due to the invalidation of the linear approximation and the beginning of a non-linear regime for density 
perturbations. The purpose of this analysis is to examine how the integrals ( p5| ) behave as the upper limit of integration 
increases. The proper inclusion of the modulating function will be performed in the next Section. 

The short-wavelength contribution to Co: In the limit of large n's, the asymptotic behaviour of Ddn) (first line in 
([lij) ) is dominated by the purely monopolar term C and the leading order expansion of term B: 

1 



£>o( 



cos(n£). 



The short-wavelength contribution to Co is given by 



Co 



47T 



100 



K(f3) I — 



1-2 



cos(n£) 



,2(/3+2) 



(35) 



This integral would be logarithmically divergent in the case f3 = —2, and power-law divergent for > —2, if one could extend 
the integration up to infinity. The short- wavelength behaviour of Co is drastically different from the one of Co, given by 
expression (p9|). 

The short-wavelength contribution to Ci: The asymptotic behaviour of D\(n) is dominated by the purely dipolar 
term A: 



£>i(r 



1 

3 n ' 



As expected, D\(n) has one extra power of n (a large number for short waves), as compared with Dd(n) and D2(n) (for D2(n) 
see formula ( |36| ) below). This makes the dipole moment much more sensitive to short waves, in comparison with the monopole 
and quadrupole, as it should. The short wavelength contribution to Ci is 
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Figure 1. The multipole distributions as described by formula (solid line) and formula ( p6| ) (dotted line) for the flat primordial 
spectrum /3 = —2 (n = 1). 



An 



,2(0+3) 



The integrand has two extra powers of n, as compared with Co and Ci- 

The short-wavelength contribution to Ct>?: The asymptotic behaviour of Dt>i[n) is dominated by term B: 



Dt>^n) 



cos(n£). 



(36) 



(Strictly speaking, this expression is valid for even £'s; for odd £'s, the cos is replaced by the sin.) The short wavelength 
contribution to Ce>2 is given by 



47T 



100 



K(f3) 



1-g 



dn 1 + cos(2n£) 2 (/3h 



(37) 



The integrand of ( |37[ ) is quite similar to the integrand of (g5|) . 

We have demonstrated the importance of the terms A, B and C in the derivation of the expected CMB multipoles. In the 
short wavelength regime, term D is always sub-dominant in Dt(n), while it is the sole term retained in formula (^) for Ct- To 
illustrate the difference between the Ct and Ct distributions, we perform an explicit numerical calculation for the flat spectrum 
(3 — —2. The integral ( poj ) is taken up to n max = 1260. It is believed that the linear approximation is valid for wavelengths up 
to 100 times shorter than the Hubble radius (see, for example, Padmanabhan ( 199S )) which gives n max ~ 100 ■ nn — 400-7T. 
The value of £ is determined by Zdcc (see equation (Q)), which we take Zdcc = 1000. 

In figure |l|, we plot the quantity £(£ + l)Ct as a function of I. (For concreteness, we take K(— 2) = 27T.) Formula ( p8| ) 
predicts a straight horizontal line (shown by the dotted line) with the monopole moment going to infinity. For the dipole 
moment I = 1, formula (12 a) gives 1(1 + i)C t \ t=1 w 1.4 • 10 4 , while formula @ predicts for i = 1 (as well as for all other ts) 
£(£+ l)Ce = 4 • 10 -2 . The ratio C\jC\ ~ 3.5 ■ 10 5 clearly illustrates the dramatic difference in the expected dipole. The higher 
order multipoles also experience a significant deviation from the straight line (mostly, due to the inclusion of term B) as shown 
by the solid line. The ratios Ct/Ct for I =2, 10, 20, 30, 40, 50 are 1.03, 1.52, 2.69, 4.35, 6.42 and 8.82, respectively. The higher 
order multipoles are more sensitive to short waves, and therefore the difference between the lines is even more prominent. The 
growth of the function £(£+ l)Ct would be unlimited, if the spectrum could be continued as a purely power-law spectrum for 
n ^> n max . However, in reality, it is terminated because of a change in the spectrum caused by the inclusion of the transfer 
function. 



5 THE PRIMORDIAL SPECTRUM AND THE MODULATING FUNCTION 

A smooth spectrum of density perturbations that existed deeply in the radiation-dominated era gets strongly modified by 
the time of decoupling. This is the result of evolution of oscillating density perturbations (metric and plasma sound waves) 
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through the transition from the radiation-dominated regime to the matter-dominated regime. This phenomenon is usually 
described with the help of the transfer function (see, for example, Peebles (1993), Efstathiou ( 199d| ), Zel'dovich & Novikov 
Q1983D , Doroshkevich & Schneider (jl99ej), Hu & Sugiyama (|l995t)) and several model transfer functions are known in the 
literature. Their common property is a strong reduction of the short-wavelength power, starting from some large critical 
wavenumber. For concreteness, we shall use the transfer function derived in Grishchuk (1994). The purpose of that paper was 
to derive the power spectrum from first principles and the quantum-mechanical approach was consistently used. In a moment, 
we shall briefly outline that derivation in order to relate notations, but the final result is not unexpected, so we shall present 
it first. The result is the following: instead of the strictly power-law spectrum 



\B„ 



K(J3) n 



2/3 + 1 



introduced by equation (E4), the integrals rt23j) for the Cis should be taken with the modified spectrum 



\B„ 



K(f3) n 2p+1 M\ni 2 ), 



where the modulating (transfer) function M (n£z) is given by 



M 2 « 2 ) = 



I "£2 \2 



(38) 



(39) 



and £2 is determined by the red-shift z cq at the transition from the radiation dominated era to the matter dominated era: 



€a = (1 + Zeq) 



-1/2 



If one does not wish to go into the details of derivation, formula (J3S|) can be treated as a postulated phenomenological transfer 
function, and formula ( |38| ) as a postulated resulting spectrum. 

We shall briefly outline the derivation of the power spectrum for density perturbations, but the reader interested only 
in applications may skip this paragraph. The theory of quantum- mechanical generation of density perturbations (Grishchuk 
1994T 



Grishchuk 1996) uses the Heisenberg operator for the field hij(r/,x): 

L {n) c n (0) e inx + hi (V) cl (0) e _ta 



hij(rj,x) = Zpi 



(2tt) 3 / 2 



d J n 



2 



Pij (n) 



(40) 



where the Planck length Ipi comes from the quantum normalisation of the field to a half of a quantum in each mode, while c n (0) 

t 3 
and cL(0) are the usual annihilation and creation operators. Two polarisation tensors Pij(i\) (s — 1,2) describe, respectively, 

a scalar and a longitudinal component of the field: 



P«(n) 



/ 2 
-iiiiij /n 



The functions h n (v) satisfy the perturbed Einstein equations and are being evolved through three successive stages of evolution: 
initial (i), radiation-dominated (e) and matter-dominated (m), with the respective scale factors 



P< -1, V < Vi, Vi < 0, 



»7c), ?7i < V < 

\2 



(41) 



a e (v) = ioa e (v 

flm(»)) = 2l H (rj~r] m y, T] 2 < T) < 77R. 
The constants participating in (^) are connected with each other by the continuous joining of a(rf) and a' (77) at the transition 

s 

points. At the (m)-stage, the growing solution for h n (v) i s 

121 2 
h„= C(n) , h n (v) = C(n) [n(r] ~ r; m )] , 

which allows one to relate this form of the solution with expressions ([]) , ([?]) . Comparing the quantum-mechanical expectation 
value of the square of the generated field hij with the ensemble averaging used in equation ([ic]), one obtains |C(n)| 2 = \B n \ 2 n, 
where the factor n comes from the extra factor l/\^2n in the definition dio|). The perturbed Einstein equations plus the 
imposed initial conditions allow one to find the spectrum |C(n)| 2 from first principles, and express this quantity in terms of 
the participating parameters and the fundamental constants (see equations (82), (81) in Grishchuk (1994) and accompanying 
discussion). The resulting function |C(n)| 2 has the form 

\C(n)\ 2 =K{p) n 20+2 A/ 2 « 2 ), 

where M 2 (n^) is given by ( |39[ ) and K(f3) is a known function of [3. For the specific case of the flat spectrum, f3 = — 2, 

K(-2) = 2n (^l) 2 . (42) 
to 
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Figure 2. The modulating function M 2 (n^2) for different values of z eq . The wavenumber corresponding to the present-day Hubble radius 
is Air. For waves longer than and comparable to the Hubble radius, the modulating function is approximately 1. At wave-numbers around 
400 the modulating function starts suppressing the primordial power spectrum significantly. 



The modulating function is an oscillatory function of n and has deep minima. This can be explained ( Grishchuk fc Sidorov 
19910as the eventual 'desqueezing' of some of the generated modes, and is also related to the standing wave pattern of the 
generated perturbations and the phenomenon of the Sakharov oscillations. 

We shall now work with the function leaving the parameters (3, K{j3), z cq and Zdoc free, to be determined by the 

comparison with observations. In figure ^, we plot the modulating function M 2 (n£ 2 ). 

In figure ^, we present the effect of the modulating function on different primordial spectra of the density contrast, 
V(k) — fc n , which in our notation is V(n) = n 2l3+5 . For sufficiently long wavelengths, the spectra retain their original form. At 
wave-numbers around 400 the spectra go through the maximum. 

The maxima of the n 2/3+5 M 2 (n£ 2 ) spectrum, that we present in figure^, appear at the wave-numbers which satisfy the 
equation 

The increase of displaces the maxima towards larger n. As an example, we present in figure ^ the displacement of the first 
maximum towards larger values of n, when one keeps z cq = 7000 and increases (3. The dotted line represents the right hand 
side of the equation (which is /^-independent), while the solid, dashed-dotted and dashed lines represent the left hand side for 
/3 = —2 , —1.9, —1.8. The position of the first maximum in the power spectrum is directly related to the position of the first 
peak in the ^-space, as will be discussed below. 

We shall now calculate numerically the multipole moments following from the formula 

47r~ 1 f°° An 2/ 



Cl = ~m K[P) / ~ M\n^)n^ + " Di(n), (43) 
J o 

which replaces formula (pj|). 



6 THE CMBR ANISOTROPY DISTRIBUTIONS CAUSED BY DENSITY PERTURBATIONS 

The modulating function makes the power spectrum oscillatory and dramatically suppresses the power at short waves n^2 3> 1, 
roughly, in proportion to n~ 2 . In figure]^, we present a typical graph of £(£ + l)Ce calculated with the modulating function. The 
chosen parameters are Zdec = 1000, z cq = 7000 and j3 = —2. The integration in formula ( fiil) is carried out up to n max = 1260 
and the quantity 2(2 + 1)6*2 (for the quadrupole, £ = 2) is normalised to 1. The dashed-dotted line shows the result without 
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Figure 3. The effect of the modulating function on different primordial spectra n 2 ^+ 5 . The dashed-dotted line refers to the original 
spectrum n 2 ' 3 "'" 5 , while the solid one to the resulting modulated spectrum n 2 ' 3 "'" 5 M 2 (n^2)- The spectra shown are: i) top left hand side, 



the fiat spectrum f) = — 2, ii) top right hand side, /3 
All the modulated spectra were calculated for z cq = 



= —1.95, iii) bottom left hand side, /3 ; 
7000. 



-1.9, iv) bottom right hand side, /3 = — 1. 



400 500 

wavenumber n 



Figure 4. The dotted line shows the function — [(ra£2)/(2%/3)] •cos((n^2)/ (2%/3)), while the solid, dashed-dotted and dashed lines show the 
function \J3 + (3/2)] ■ sin((n^2)/(2\/3)) for /3 = —2, —1.9, —1.8 respectively. The red-shift at matter-radiation equality is taken z CCL = 7000. 
The values of the wavenumber where the intersection of the curves takes place (except the point n = 0) are the points where the 
modulated power spectrum n 2 0+5 M 2 (n^2) has its first maximum. For increasing /3, the first maximum occurs at larger values of n. 



including the modulating function. The solid line clearly shows the first peak formed due to the inclusion of the modulating 
function. 

Before proceeding to other choices of the parameters, it is instructive to explore the integrand n 2l3+3 M 2 (n$2)D{(n) in 
formula (p3[). We shall do this for various £'s and at different intervals of integration over n. We shall present the graphs of 
the integrand, which allow one to evaluate visually the surface area under the graphs and, hence, the contributions to the 
integrals accumulated at different intervals of n. 

I = 0: In figure [], we plot the integrand of Co at different intervals of n. A deep minimum of the integrand caused by the 
first zero of the modulating function is clearly seen. It is obvious from the graph that the integral is mostly accumulated at 
small n's, up to n ~ 12 or so. Thus, the dominant contribution to Co comes from wavelengths larger than and comparable to 
the Hubble radius. 

I = 1: In figure 0, we plot the integrand of C\. The modulating function suppresses the short- wavelength power, but still 
leaves the dipole moment Ci much more sensitive, than Co and C2, to short waves. The remarkable fact is that the main 
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contribution comes from wavelengths about 20-40 times shorter than the Hubble radius. For the conventional values of the 
Hubble parameter, this corresponds to scales around (100 - 200) Mpc, not much larger and not much smaller. The most 
important term in the decomposition of Di(n) (see formula ( |l4] ) ), is term A. It is obvious from the graph that the numerical 
value of the dipole is much greater than the value of the monopole and (as we shall show below) the £ > 2 multipoles. 

I > 2: There exists a qualitative difference between small multipoles i « 2—10 and relatively large multipoles £ ~ 100—200, 
in the sense of the different roles played by terms B and D (see the decomposition of D^^n) in formula ( |l4| ) ). In the case of 
small multipoles, the dominant contribution to the integral comes from long and intermediate wavelengths, where the most 
important term in De>2{n) is term D. In the case of relatively large multipoles, the dominant contribution comes from short 
wavelengths, where the most important term is term B. In figures g, [] and [n], we present by a solid line the behaviour of the 
total integrands for C2, C10 and C200 at different n's. The values of the parameters are z cq = 7000, Zdcc = 1000 and j3 = —2. 
At the same graphs, we have also plotted with a dashed-dotted line the integrand that would result from neglecting term B, 
and with a dotted line the integrand that would result from neglecting term D. 
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Figure 7. The integrand of the dipole moment 1=1 versus the wavenumber n. The values of the parameters are again (3 = — 2, 
z cq = 7000, z dcc = 1000. 

From figures ^| and [], it is evident that the dominant contribution to Ci and Cio comes from wavelengths longer than 
and comparable to the Hubble radius. In this regime, the role of term B is unimportant. Although term B will eventually 
dominate in the short wave regime, the contribution of that regime is negligible. At the same time, the modulating function 
in the most sensitive part of integration is close to 1. Therefore, the inaccurate formula ( pj| ) gives roughly correct numerical 
values for small multipoles. On the other hand, as figure [To] shows, the main contribution to the relatively large multipoles 
comes from short scales, where term B is always dominant. The neglect of term B gives wrong results. 

Summarising, one can say the following. The major contribution to the monopole moment comes from wavelengths larger 
than and comparable to the Hubble radius. The actually measured isotropic temperature To contains this small contribution. 
The major contribution to the dipole moment comes from wavelengths significantly shorter than the Hubble radius, but still 
far away from the onset of nonlinearities. Numerically, the dipole moment is much larger than the monopole and the I > 2 
multipoles, as seen from the comparison of figures ^ and p|-|lo| with [?]. To illustrate the greater sensitivity of the dipole to the 
short waves, one can consider the specific example of (5 = —2, z cq = 7000 and Zdcc = 1000. The ratio of the dipole to the 
quadrupole is C1/C2 ~ 1.4 • 10 5 . This ratio in the absence of the modulating function is ~ 10 6 . The major contribution to 
the t € [2, 20] multipoles comes from wavelengths longer than and comparable to the Hubble radius. The major contribution 
to the relatively large multipoles comes from short wavelengths, where term B dominates. Still, these multipoles are several 
orders of magnitude smaller than the dipole. 

6.1 Variation of the parameters 

We have demonstrated which wavelength regimes are responsible for the main contribution to different multipoles. While 
calculating the integral (^), one particular (typical) choice of the parameters has been used: Zd ec = 1000, z cq = 7000 and 
P = —2. One needs to explore how the change of these parameters affects the CMBR anisotropy pattern. The factor Di(n) 
depends on Zdec through the presence of £. The modulating function depends on z eq through the presence of £2- The primordial 
spectrum depends on /3. We shall separately discuss the effect of the variation of Zdec, z cq and 0. The results are presented in 
figure |ll| 

The effect of Zd ec : The values of the fixed parameters are z eq = 7000 and /3 = —2, while we decrease Zdec from 1000 
(solid line) to 500 (dashed-dotted line). The graph shows that the peak is displaced towards lower multipoles £, while the 
height of the peak increases. As already stated, the decrease of Zdec decreases £. This leads to the increase of term B, and 
explains the rise of the height of the peak. The dipole moment is insensitive to the change of Zdcc, since it is dominated by 
term A which is not affected by Zdec- 

The effect of z eq : The values of the fixed parameters are Zdec = 1000 and j3 = —2, while we allow z eq to be 5000 
(dashed-dotted line). The graph shows that the peak is displaced towards lower multipoles £, while the height of the peak 
decreases. The decrease of z oq increases £2, making the modulating function decline faster (figure |^). The Ce will acquire 
smaller values, and therefore the values of 1(1 + l)Ce will be smaller. This accounts for the decrease of the height of the 
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Figure 8. This plot shows by a solid line the integrand for the quadrupole moment, I = 2, at different wavenumber intervals. We also 
display the integrand that results after retaining only term D in D^in) (dashed-dotted line), and the integrand that results after retaining 
only term B (dotted line). For long wavelengths, term D clearly dominates over term B, while for short wavelengths the situation is 
reversed. At n ~ 900, we observe a zero caused by the first zero in the modulating function. The modulating function and the Bessel 
functions combine in such a manner, that the area under the curve at long and intermediate wavelengths is much greater than the area 
at short wavelengths, where term B dominates. The main contribution to C2 is provided by large scales. 
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Figure 9. The same as in figure |], for the multipole I = 10. Although the role of term B is now more important than in the case of the 
quadrupole, the main conclusions of the i = 2 case still apply to I = 10. 
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Figure 10. The same as in figure g, for the multipole t = 200. The main contribution to this multipole comes from relatively short 
waves, where term B dominates. The modulating term and the Bessel functions combine in such a manner, that the area under the 
curve at long and intermediate wavelengths is much smaller than at short wavelengths. The situation is opposite to the case of small 
multipolcs. 



peak. At the same time, the decrease of z eq displaces the first minimum of the modulating function (and of the modulated 
spectrum) to smaller n's (figure ^). The first maximum of the modulated power spectrum will then occur at a lower value of 
n, which translates in the displacement of the maximum of £(£ + l)Ct towards smaller £'s. This is the dominant reason for the 
first peak displacement, but we do not exclude that some additional effects can counterbalance this tendency. The decrease 
of z eq decreases the value of the dipole. The ratio Ci|z,. q =7ooo / Ci|z =5000 = 1.25. 

The effect of j3: The values of the fixed parameters are z cq = 7000 and Zdec = 1000, while we allow /3 to be —1.9 
(dashed-dotted line). The graph shows that the peak is displaced towards higher multipoles I, while the height of the peak 
increases. The displacement of the peak is explained as follows. As shown in figure ^, the increase of f3 displaces the position 
of the first maximum of the modulated spectrum towards larger n. As already discussed, the change of f3 does not affect the 
position of the minima. Therefore, the maximum of the function £{£ + l)Ct versus £ is displaced towards higher t. As for the 
increase of the height of the peak, it is explained as follows. The primordial flat spectrum n 2/3+5 (top left graph of figure ^, 
dashed-dotted line) is a smooth line, growing with n. For a larger fl, the primordial spectrum increases faster with n, and 
its plot versus n is a concave curve. Therefore, the amplitudes of the modulated spectrum are larger for a larger /3, which 
increases the values of the Ce and the height of the peak. The increase of (3 increases also the value of the dipole. The ratio 
Ci\g=-i.9 I Ci\g = -2 « 3. 



6.2 The inclusion of intrinsic anisotropies 



The derivation of the starting formula ^ is based on solving the differential equation for the photons' geodesies (Sachs & 



Wolfe 1967). As usual, when defining a solution, one should indicate the initial conditions. In our problem, this translates 
into the intrinsic CMBR anisotropies. These anisotropies reflect the temperature inhomogeneities already present on the last 
scattering surface. There are many elaborate treatments dealing with this subject (for example, Doroshkevich & Schneider 
(1996), Hu & Sugiyama ( |1995| )). We shall use a simplified model (see, e.g. Padmanabhan (11993)), where the intrinsic anisotropy 
is related to the radiation density contrast by the simple relationship 
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Figure 11. The effect of change of different parameters. The quantity 1(1 + l)Ci is normalised to 1 at £ = 2. In all graphs, the solid 
line corresponds to z^ cc = 1000, z eq = 7000, (i = —2. i) Top graph: the effect of decreasing z<jcc- The dashed-dotted line corresponds 
to 2(j cc = 500, z cq = 7000, P = —2. The peak is displaced to the left and its height increases, ii) Middle graph: the effect of decreasing 
z cq . The dashed-dotted line corresponds to z^ cc = 1000, z cq = 5000, j3 = —2. The peak is displaced to the left and its height decreases. 



iii) Bottom graph: the effect of increasing (3. The dashed-dotted line corresponds to z^ e . 
displaced to the right and its height increases. 
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Deeply in the radiation-dominated era, photons are tightly coupled to baryons through Thomson scattering. It can be seen 
that the radiation density contrast can be related to the baryon density contrast by 



Assuming that the dominant matter component governing the scale factor at the matter-dominated stage is a sort of dark 
matter, the baryon density contrast can be related to the dark matter density contrast by 



-ec 



Therefore, one can relate the intrinsic anisotropics to the total density variation at the last scattering surface: 



(f), 



1 -Zdc 
3 -Sec 



This intrinsic part should be added to the gravitational part of the previous Sections, and this total temperature variation 
should be decomposed over spherical harmonics. One then repeats the procedure already performed in the absence of the 
intrinsic anisotropy, to derive the full formula (containing both gravitational and intrinsic parts) for the multipole moments 
Ci. This is again equation (fl3|), where the expression ( |l3| ) for Di(n) is modified to 



Di(n) = - n Sn + n(l 







■[(£+1) 3t+i(nQ-tjt-M)] + 8tD 



1 + — - (nf 2 ) 
10 6 v s ; 



(44) 



2£ + l' 

The inclusion of the intrinsic part causes the further increase of the multipole moments Ce, and the increase of the height of 
the first peak in the CMBR distribution. If plasma fluid at the last scattering surface retains appreciable velocity with respect 
to chosen by us coordinate system, which is synchronous and comoving with gravitationally dominant matter, this will give 
rise to further corrections. As will be argued in the next Section, the observed first peak seems already to be too high to be 
explained, together with other multipoles, by density perturbations alone. 
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7 A FIRST COMPARISON WITH OBSERVATIONS 

The current CMBR observational results can be summarised as follows. 



1) The COBE observations ( Bennett et al. 1996] ) have shown the presence of a 'plateau' in the multipole region £ G [2, 20]. 
As a characteristic value in this region, we shall use the quantity £{l + l)C^|<=io. By the bold-face symbol Ce we denote the 



value of the multipole moment £ that follows from observations. (The lack of ergodicity on a 2-sphere ( Grishchuk & Martin 
199?|)]prevents one from precise extraction of statistical Ci from a single, even if arbitrarily accurate, map of the sky, but this 
is not a matter of concern for us here.) The COBE result shows that 



t(! + i)c t \t=i 



271 



To = 27.7±1;1 fiK 



where To is the background, isotropic, temperature of the CMBR. [|] This gives the value of £{£ + 1)C^ 



l)C*|i=xo « 6.4- 10" 



(45) 



2) The dipole anisotropy was also measured by COBE. According to Bennett et al. (1996) 



(2l + l)C £ |i=i 

leading to 

Ci w 6.3- 10" 6 . 



To = 3.353 ± 0.024 mK, 



(46) 



3) All the CMBR observations performed at the region of higher £'s have suggested the presence of a peak. The most recent 
and accurate experiments BOOMERanG (de Bernardis et al. 2000) and MAXIMA-1 (Hanany et al. 200C) gave the location 
and height of this peak with some precision. According to these results, the peak is located at 200 < ^ pca k < 220, with 
y/£{i+ l)C f | pcak /(27r) • T having the value 69 ± 8 fM ( |de Bernardis et al. 2000| ) or 78 ± 6 fiK ( [Hanany et al. 2000[ ). For the 
sake of this discussion, we adopted the values of Hanany et al. (200C), 



^peak 

and 



220 



(47) 



l(l+l)C,| pcak 
2tt 



To = 78 ± 6 fj,K. 



Therefore, we shall use as a characteristic value in this region the quantity 
l(£+ l)Cf| peak « 5.1 -10" 9 . 



(48) 



We accept the position of the peak (^), and we regard the three numbers (jig), (pf^), ( f45| ) as three observational points 
to which we try to fit (admittedly, in a somewhat simplistic manner) the theoretical curves. The curves for £{£ + l)Ce are 
constructed with the help of formulae ([|3]) and (^), where the upper limit of integration is taken at n max = 1260. The 
strategy of the fitting is as follows. 

There are several combinations of the parameters (.s C q, z^cc-, (3) (some of them are listed in table |l|) for which the position 
of the peak is at © and the height of the peak is at (^). At the same time, we check that the calculated statistical dipole is 
close (deviation less than 5 per cent) to the observed value (|46|). A drastic adjustment of the calculated dipole with the help of 
a 'peculiar' velocity of the individual observer, that is, allowing the observer to move with a large velocity with respect to the 
matter-dominated fluid (i.e. with respect to the comoving coordinate system) seems to be artificial. As the statistical dipole is 
already quite large, it seems especially unnatural to try to do this. As a guideline in the choice of Zd ec an d 2 eq , we use the fact 
that 2dcc should be close to 1000 and -z eq should be somewhere near to 10 4 . The decoupling red-shift Zdcc certainly depends on 
the baryon content of the Universe fib, but it can also depend on a complicated ionisation history of the decoupling era. The 
equality red-shift z eq certainly depends on the total matter content of the Universe (baryons, cold dark matter, etc.) but it 
also depends on its radiation content (CMB photons, possible contribution of massless neutrinos, etc.). In addition, the total 
density parameter Q can also deviate from the postulated here £7 = 1. This may change z cq without seriously affecting our 
solutions for the gravitational field perturbations hij and subsequent formulae, derived, strictly speaking, for the case Q = 1. 
It would be a difficult problem to try to decipher all the relevant quantities from just the two numbers z^ cc and z cq . However, 
we will operate with these parameters, since this is what the calculations immediately require us to specify. As for /3, we have 



t In many topics related to CMBR. experime nts and observational data, we have resorted to the web-site 



http: / / www. hep. upcnn.cdu~max/cmb/cxpcriments. html. 
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Figure 12. Comparison of theoretical models, based on density perturbations only, with observations. The solid line corresponds to the 
set of parameters (5000, 1000, -2), the dashed-dotted to (4500, 1000, -1.95) and the dotted to (4000, 1000, -1.9). The crosses show the 
observed values of t{i + l)C|^ = io and £(£ + ^C^p,,^. By fitting the theoretical peak to the observed peak, we notice that there is a 
deficit at the plateau. 



Table 1. Comparison of some theoretical models, strictly satisfying the requirements on the position and the height of the first peak, 
with observations. 



•2-eq 


z dcc 






A 




£{£ + l)C e \ e=10 


Cio/Cio 


Ci 




Ci/C 


5000 


1000 


-2.00 


7.7 


10~ 


-9 


4.2 ■ 10" 10 


1.5 


6.1 ■ 10" 


<i 


1.04 


4500 


1000 


-1.95 


5.0 


10~ 


-9 


3.4 ■ 10" 10 


1.9 


6.1 ■ 10~ 


-6 


1.03 


4000 


1000 


-1.90 


3.3 


10~ 


-9 


2.9 ■ 10 -10 


2.2 


6.1 ■ 10~ 


-6 


1.03 



considered the three different cases f3 = —2, —1.95, —1.9. Since K(/3) is not very sensitive to (3 in this interval, we have used 
one and the same formula ([42]) for all the three cases. Then, every fit of the peak and the dipole to their observed values defines 
its own value for the ratio A = (Zpi/io) 2 and for the calculated quantity £(£ + l)Ce\e=io- This quantity is finally compared 
with the observed value (j^f. The result of this comparison shows that the theoretical £(£+ l)Ce\e = io is systematically smaller 
than the observed £(£ + l)C|f = io, for all considered combinations of the parameters (z cq , Zdcc, /3). The observational points 
should be surrounded by their error bars, but in order not to proliferate the already existing uncertainties, we focus on the 
mean values. One can say that, despite of the still allowed variation of the parameters, the calculated 'plateau' is persistently 
lower than it should be. These calculations are visualised in figure 

We believe that this deficit at small multipoles serves as an indication that density perturbations alone can hardly be 
responsible for all of the CMBR anisotropies. The most natural way to proceed is to include gravitational waves and repeat 
this scheme. 
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8 THE GRAVITATIONAL WAVE CONTRIBUTION TO THE CMBR ANISOTROPIES 



Gravitational waves are inevitably generated by the strong variable gravitational field of the early Universe (see, for example, 
Grishchuk et al. ( 2000) and references therein). The preliminary calculations show (Grishchuk 1994) that the contribution 
of quantum-mechanically generated gravitational waves to the small multipoles should be of the same order of magnitude 
and, numerically, somewhat larger than the contribution of density perturbations generated by the same mechanism. In this 
paper, however, we follow a phenomenological approach and do not not focus on the origin of cosmological perturbations 
and consequences of the particular generating mechanism. We introduce, instead, a phenomenological parameter zu, which 
regulates the relative contributions of gravitational waves and density perturbations. In effect, we want to find the value of vj 
which allows to raise the plateau in figure jl2| to the observationally required level, at the expense of gravitational waves. At 
the same time, the fits to the peak and the dipole are required to be as good as before. It is already clear from table |l| and 
figure [l2] that the required amount of gravitational waves, as compared with density perturbations, is within a factor of 2 at 
£ — 10. This seems to be in agreement with theoretical expectations, but we shall go into a more detailed analysis in Section 



In this paragraph, we shall only remind the fundamentals about primordial gravitational waves. Although we follow the 
quantum mechanical derivation, the final formula (^lj) for Ci is true (up to the overall numerical coefficient) for a broad class 
of gravitational waves. The Heisenberg operator for the gravitational wave metric perturbations is 



hij(ri,x) = Zpi 



(2tt) 3 / 2 



2 

E 

3=1 



Py(n) 



h n {v) c n (0) e in x + h* n {n) 4(0) e- ir 



(49) 



where, again, Zpi comes from quantum normalisation, and the transverse traceless polarisation tensors (n) satisfy the 
relations 



Pij(n) p ij (n) = 2S SS ,, Pij{n) =p ij (-n), P i3 n j = 0, P tj S lj = 0. 

The time evolution of each of the amplitudes h n (rj) for s = 1, 2 is given by one and the same equation with the same initial 
conditions, so one deals with one and the same h n (rj). After solving for h n (r)) at the three stages ( |4lj ) and continuously joining 
the solutions, one arrives at a 'growing' and a 'decaying' solution in the matter-dominated era. The growing solution is given 
by the expression 



ftn(i))wie in "» 3*(/3) (-1)- 



,-1 3+1 h(n(v -Vm)) 



(50) 



where the quantity ^(/J)! 2 is unity for /3 = —2, and we neglect its weak /3-dependence for other values of (3 considered here. 
For a reference, [\&(/3)| 2 is related with K(j3) discussed in Section ^ by 

m = **\m\ a (jff- 

Combining (^o|) with (^), using the Sachs- Wolfe for mula ([jj), and re peating the same steps that have led to the formula for 
Ce's caused by density perturbations, one arrives at (Srishchuk 1993) 



d ^4-9|*(/3)| 2 
where 
h{n) = 



jix) jdn -x) 



(n — x) 



l)t{t +!)(£ + 2) 



dx 



2(3+2) r I \ 

^ yp ' fe(n) — , 



(51) 



As ( |5l] ) shows, gravitational waves contribute to multipoles £ > 2, in contrast to primordial density perturbations, which 
contribute to all multipoles. 

We shall now present the multipole moments caused by gravitational waves. The solid line in both graphs of figure |l3"[ 
shows the behaviour of £(£ + l)Ci versus £ for z cq = 7000, z^ cc — 1000 and j3 — —2. The quantity £(£ + l)Ce\e=2 is normalised to 
1. The higher order multipoles are negligible in comparison with the lower order ones, unlike the case of density perturbations. 
This is expected because of the different behaviour of the short wavelength density perturbations and gravitational waves. 
At the matter-dominated stage, short wavelength density perturbations increase with time as rf , while short wavelength 
gravitational waves decrease with time as r\~ 2 . Figure [l3] also shows the effects of the variation of different parameters. As 
expected, in all cases, gravitational waves mostly contribute to the small multipoles. 
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Figure 13. The effect of the variation of the different parameters on 1(1 + l)Ci. In both graphs, the solid line corresponds to Zdcc = 
1000, z C q = 7000, f3 = —2. i) Left graph: the effect of decreasing Zdec The dashed-dotted line corresponds to z^ cc = 500, z cq = 
7000, /3 = —2. ii) Right graph: the effect of increasing /3. The dashed-dotted line corresponds to z^oc = 1000, z cq = 7000, /9 = —1.9. 
The local maximum is displaced to the right and its height increases, exceeding the value of £(£ + l)Cf \i—2- 



Table 2. Comparison of the theoretical model with observations, for the sets of parameters (5000,1000,-2), (4500,1000,-1.95) and 
(4000, 1000, —1.9). Primordial gravitational waves have now been taken into account. The position and the height of the first peak, as 
well as the dipole moment and the plateau, are in an excellent fit to the observations. 
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3.4 ■ IO" 10 
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9 A SECOND COMPARISON WITH OBSERVATIONS - CONCLUSIONS 

We denote by Cf v and C| w the multipole moments induced by density perturbations and gravitational waves, respectively. 
Cf p is given by (^), whereas Cf w is given by (|5ll). We introduce a numerical parameter w in such a way that the total 
£(£+l)Ct is given by 

lit + l)C t = £{l + l)[Cf w + zu ■ C^]. (52) 

The value of vj would be strictly 1, if one could guarantee that the idealised model that has led to formulae (fl3|), has 
accurately taken into account all physical processes that might have affected density perturbations and gravitational waves 
throughout their entire evolution, from the moment they existed in a vacuum quantum state and up to the era of decoupling 
of CMB radiation from matter. Since one may have legitimate doubts in such accuracy, the correction parameter w may be 
thought of as a measure of a necessary modification (but only in the form of an overall numerical coefficient) to our theoretical 
framework. 

As was shown in Section gravitational waves do not contribute to the dipole, and do not practically contribute to the 
peak. Therefore, we require the quantity zu ■ £(£+ l)Cy P j pca k to be equal to the characteristic value (^]). At the same time, we 
require £(£ + l)CVjf=io, defined by setting £ — 10 in (p^), to be equal to the characteristic value (^). This comparison revises 
the ratio A = (ipi/io) 2 and determines the parameter vj. We then construct the ratio £(£+ l)Cf w / (zu ■£(£+ l)C^ p |£=io) = 
Cfo" / (zu-Ciq), which shows the relative contributions of gravitational waves and density perturbations to the small multipoles. 
Finally, we compare the value m ■ Ci with the observed value ( [l6|) of the dipole. The results of this analysis are shown in table 
|j| while a visual representation is given in figure 
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I 



Figure 14. Comparison of the theoretical model, including density perturbations and gravitational waves, with observations. The dotted 
line shows the contribution of the corrected density perturbations alone, the dashed-dotted line shows the contribution of gravitational 
waves alone, and the solid line is the total £(£ + 1)C£. The top graph corresponds to the set of parameters (5000, 1000, —2), the middle 
graph to (4500,1000,-1.95), and the bottom graph to (4000,1000,-1.9). The crosses, again, correspond to the observed values of 
1(1 + l)C|f— io and 1(1 + l)C| pcak . In all cases, the contribution of gravitational waves to the plateau is comparable to that of density 
perturbations, and for larger values of /3 the contribution of gravitational waves is numerically greater than that of density perturbations. 



To summarise, we have revised the formula (^) for Cts caused by density perturbations and derived the corrected 
expression (fl3|). We have shown that the monopole moment Co is actually finite and small, while the dipole moment is very 
sensitive to short wavelengths. Because of that, the dipole moment is much larger than the I = and I > 2 multipoles. We 
have then proceeded to the multipoles I > 2. The correct formula predicts growth of the function til + l)Ce with t. We 
have illustrated how the presence of the modulating function M\n^2) gives rise to the first peak in the CMBR anisotropy 
pattern produced by density perturbations. (It would be nice to see further peaks discovered in the positions and with the 
heights suggested by the modulating function, but this is the domain of wavelengths where many other physical processes 
may intervene.) A first comparison with observations has indicated that the observed multipoles can hardly be explained by 
density perturbations alone, because of the deficit in the small £ plateau. We took a step further and included primordial 
gravitational waves. This inclusion gives a picture in better agreement with observations. We explored the cases > —2 which 
make the primordial spectrum a little bit 'blue', thus allowing to avoid the infrared divergence of the primordial perturbations 
that occurs for < -2 (see @ and @). We conclude that models with > —2, and therefore with larger contributions of 
gravitational waves, are better than a model without gravitational waves and with the flat spectrum = —2. This analysis is 



complementary to other studies of the CMBR data (Albrecht 200C; Balbi et al. 2000a; Balbi et al. 2000b; Durrer & Novosyadlyj 


200C; 


Hu ct al. 2000; Kinney 200C; 


Lange et al. 2000; 


VtcGaugh 200C ; Martin, Riazuelo & Schwarz 200C ; Padmanabhan & 


Sethi 


2000; Peebles, Seager & Hu 2000; Pogosian 200C 


; Rahman & Shandarin 2000; Schmalzing, Sommer-Larsen & Goetz 


200C; 


Tegmark & Zaldarriaga 200C; 


Tegmark, Zaldarriaga & Hamilton 2000a|; Tegmark, Zaldarriaga & Hamilton 2000b), not 



all of which, though, seem to agree with each other. 
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APPENDIX A: OUTLINE OF THE CALCULATION OF THE CONTRIBUTION OF TERM II 

We need to prove that the contribution of term II to formula (b|) is given by 



4tt 



(2tt)3/2 



d 3 n Bn n(l - £) 



2£ + l 



[(e + 1) - 1 ji-M)] Yfje, 



(Al) 



The starting point of the calculation is the expression 
1 f +oc 

Gll{e) = (2^7* J x ^ & H) " ' 6 (1 * 
which needs to be presented as 

Gn(e) = J2 E g£m Y ^,4>)- 

t=0 m=-£ 

Using the representations (^) and (|l^) in the scalar product n • e, the spherical wave expansion of a plane wave ( |rij ) , and the 
explicit form of Yif(9, <f>) (see e.g. page 99 of Jackson ( 1975)) we derive: 



(A2) 



(A3) 



=0 m—-t /=-! 



in • e (1 - e ine « = 4tv ^ ^ ^ i t+1 ^ n(l - £) YJ> (0, *) F/ m (G, #) F x/ (0, 0) Y im (d, 



(A4) 



We now equate the right hand sides of (A2) and (A3), multiply the resulting expres sion with Yp m r(6,tj>), and use the 
orthogonality condition of the spherical harmonics (e.g. formula (3.55) in Jackson ( 1975] )). This allows one to write 



12 3 
W = &l'm' + Sl'm' + &l'm' 



(A5) 



where 



4ti 



e=o m=-l 



d 3 n Bn n(l - #n0 ^-1(6, S) F/ m (0, $) 



Yi,_i(0, (/.) Y^(M) Y^* m '(M) sin0 dO d0 



(A6) 



4-7T 



47T 



(27r) 3 / 2 ^ ^ 3 
£=0 m=-« 



EE 



d 3 n Bn n(l - jX<) n*o(e, $) F/ m (G, *) 



Yi.de, </>) Y lm (6, <j>) Y e *, m ,(8, <f>) sinO d6 d<f> 



(A7) 



4ti 



(27T) 



3/2 



EE 



4?T .£ +3 

3 



d 3 n a, n (i - o j,(<) n*i(e, $) r; m (e, $) 



Y m (0,<« y/w(M) sine d6» <ty. 



(A8) 



In each of (A6)-(A8), we decompose the spherical harmonics into associated Legendre functions (formula (3.53) of Jackson ( 



1975)), perform the integration over 6 and <f>, utilising formulae (8.733.2)-(8.733.4) of Gradshteyn & Ryzhik (198C) and the 



orthogonality condition of the associated Legendre functions (formula (3.52) of Jackson ( 1975)), and find these expressions 
to give 



47T 



■e'+2 



( 2vr )3/2 
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a/1 — B 2 n m'+l/t2\ i • / y-\ Vl — B 2 n m' + tfr}\ 

o p e>+i ( B ) + Jt'-iiO o P i'-i ( S ) 



(A9) 



4% .£> +2 



(2vr)3/ 



d n Bn n(l — ?J 



4?r (£' + m')! 2^' + l 



j//+i(n€) [-(£' -m' + 1)] S P^ +1 (B) + jv-M) V + m') B P^(B) 



im'<E> 



(A10) 



4tt .i' +2 



(2tt)3/ 2 
J>+l(«£) 

where we have used the shorthand notations 
A = cos#, B = cosO. 



d 3 n Bn n(l-£) ' 



4tt (£'+m')! 2£' + l 
C^-"' + !)(<' -"' + 2) + 



-(£' + m' - l)(f + m') 



\/l-B 2 P™'_-\B) 



-im' <& 



(All) 



Substituting (A.9)-(A11) in (A5), and using formulae (8.733. 2)-(8.733.4) of Gradshteyn & Ryzhik (1980) to show that 



P%i\B) — (t — m + 1) B P? +1 {B) -l(e- m + l)(£ -m + 2) ^J\-B 2 P^\B) = -{I + 1) Pr(B) 



and 



lVl-B 2 Pri\B) + (£ + m)B Pr-i{B) -l(£ + m- 1)(£ + m) \/l-£ 2 P^\B) = £ PP(B), 



we finally arrive at (Al) 
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